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ABSTRACT 

We measure the binarity of detached M-dwarfs in the Kepler field with orbital periods in the range 
of 1-90 days. Kepler's photometric precision and nearly continuous monitoring of stellar targets over 
time baselines ranging from 3 months to 4 years make its detection efficiency for eclipsing binaries 
nearly complete over this period range and for all radius ratios. Our investigation employs a statistical 
framework akin to that used for inferring planetary occurrence rates from planetary transits. The 
obvious simplification is that eclipsing binaries have a vastly improved detection efficiency that is 
limited chiefly by their geometric probabilities to eclipse. For the M-dwarf sample observed by the 
Kepler Mission, the fractional incidence of eclipsing binaries implies that there are O.lllg Q 4 close 
stellar companions per apparently single M-dwarf. Our measured binarity is higher than previous 
inferences of the occurrence rate of close binaries via radial velocity techniques, at roughly the 2a 
level. This study represents the first use of eclipsing binary detections from a high quality transiting 
planet mission to infer binary statistics. Application of this statistical framework to the eclipsing 
binaries discovered by future transit surveys will establish better constraints on short-period M+M 
binary rate, as well as binarity measurements for stars of other spectral types. 

Keywords: stars: binaries: eclipsing - stars: binaries: close - stars: low-mass - methods: obser¬ 
vational - methods: statistical - techniques: photometric 


1. INTRODUCTION 

Our existence in a planetary system around an isolated 
star makes it easy to forget that single, G-type stellar 
systems are not the norm throughout the Galaxy. For 
example, among the FGK dwarfs in the Solar nei ghbor- 
hood, only 54% are single dRaghavan et al.1 [201011 . Fur¬ 
thermore, studies of the s tellar populations of open clus¬ 
ters (e.g. iGhabrierl 120031) and the Solar neighborhood 
(jHenrv et al.1 120061 iWinters et all 120151 ) have demon¬ 
strated that over 70% of the hydrogen-burning stars in 
the Galaxy are in fact M-dwarfs, i.e. stars whose masses 
are below O.GAfg. Multiple-star systems and low-mass, 
M-type stars are therefore common outcomes of star for¬ 
mation. 

Observed metrics of multiplicity across various stel¬ 
lar types provide important constraints for star forma- 
tion and dynamical ev olution theories (see review by 
iDuchene fe Kraual2013l and references therein). For the 
burgeoning field of exoplanets, however, there are several 
other reasons to study the binarity of M-dwarfs. 

For one, it turns out binary stars can host plane¬ 
tary systems, a Cen B is a particularly noteworthy ex¬ 
ample, with an Earth-like planet candidate in a short- 
period orbit around a member of a triple-s tar system 
i|Dumusaue et al.ll2012t [Plavchan et al.l 120151) . In fact, 
it may be quite common for planets to orbit one com- 
ponent of a wid ely-spaced binary system. For instance, 
iNgo et all 1 2015 ) find that, for stars hosting hot Jupiters, 
the incidence of stellar companions is ~ 50%. Planets 
have also been discovered in circumbinary configurations. 
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At first considered only exotic possibilities, circumbinary 
planets are a new class of exoplanetary system discovered 
by the Kepler Mission, wi th one or more planets orbiting 


exter i or to a pair of stars dDovle et al 1 120 111 


2012 ; lOrosz et al.l [2012 1. Rostov et al.l 1201 


Welsh et al 
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2014). 


M-dwarfs are also prolific planet hosts, particularly 
for planets with radii and masse s less than those of 
Neptune. iMorton fe Swift! (|2014fl find that there are 
on average 2.00 ± 0.45 planets per M-dwarf with peri¬ 
ods less than 150 days, based on the statistics of Ke¬ 
pler transit detections, with a peak in the radius dis - 
tribution near ~ 1 R®. iDressing fe Charbonneaul (12015D 
find that, conservatively, there exist 0.18 Earth-sized 
planets per M-dwarf habitable zone. Owing to their 
favourable brightness contrast and the close-in loca¬ 
tion of their habitable zones, planets around M-dwarfs 
offer great potential for discovery and characteriza¬ 
tion with current and f utur e miss ions, such as MEarth 
dNutzman fe CharbonneaiTl 120081) . K2 (|Howell et al.1 
120141) . TESS dRicker et ah 2014), and u pcom ing RV sur¬ 
veys (cf also IDressing fc Gharbonneaull2013 ). 

Binaries can also confound the effort to search for plan¬ 
ets. An eclipsing binary system in the background of a 
target star could mimic foreground planet transit signals, 
forming an important source of astrop hysica l false pos¬ 
itives for transiting plan et surveys (|Torres et al.l 120111 : 
IMorton fc Jolmsonll201ll ). 

Given these considerations, M-dwarf multiples con¬ 
stitute an important subclass of objects for stellar 
and exoplanetary science. The first large multiplicity 
study targeting low-ma ss primaries was conducted by 
iFischer &; Marcvl (I1992D (hereafter FM92). In this sem¬ 
inal investigation, the authors combine surveys based 
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on a variety of techniques, including radial velocities, 
high-contrast imaging, and interferometry, each sensi¬ 
tive to a particular orbital separation regime. In their 
calculation of multiplicity statistics of systems with M- 
dwarf primaries, FM92 measure frequency distributions 
in semimajor axis and mass ratio. Overall, they find that 
42%±9% of the M-dwarfs studied are in multiple systems 
within the separation range of 0.04 to 10 4 AU. 

Since FM92, the occurrence statistics of low-mass stel¬ 
lar systems have been the subject of various works. 
Some highlights are: orbital elements of spectro- 
scopic binaries among local early- to mid-M-dwarfs 
(|Mavor et al.1 120001 : lUdrv et all I2000D, volume-l imited 
study of M-dwarfs within 9pc dDelfosse et al.ll2004f) . latc- 
M-dwarf primaries (lAllenll2007tl. visual M-dwarf systems 
dBerg fors ct al ■ 10.10; .T anson c t al . 120121 I.Todar et~abl 
1201, it iWard-Duong et al.1 I2015D . young M-dwarf pri¬ 
maries ( Bowler et al 7 120151), X-ray bright M- dwarfs and 
spectroscopic binaries dShkolnik et al.ll2010D. and very 
short -period M-dwarf binaries dClark. Blake, fc Knapd 
120121) . These studies rely on either direct imaging or 
spectroscopic/radial velocity (RV) techniques. The for¬ 
mer has limited sensitivity to low projected component 
separations and contrasts, whereas the latter theoreti¬ 
cally makes up for completeness for very closely orbiting 
systems, but requires bright targets. The overall multi- 
plcity for M-dwarf primaries found by these studies lie 
between 20% and 40%. 

The closest separation bin considered in FM92 is 
0.04 AU to 0.4 AU. For the typical M-dwarf masses in¬ 
volved, such a range of semimajor axis corresponds to 
orbital periods of approximately P < 90 days. This 
very tight orbital range of M-dwarf binaries is interest¬ 
ing yet enigmatic territory. It is interesting in provid¬ 
ing constraints to theor i es of tight binary formation (e.g . 
IBonnell fc Batel Il994bl : IBate. Bonnell fc Bromml 120021) . 
It is also a discovery space for transiting circumbinary 
planets, for which the binary orbital period needs to lie 
well within the orbital period of any detectable planet. 
Currently, the sensitivity for the latter is no more than 
a couple of years, limited by survey duration. Further¬ 
more, background EBs with periods comparable to the 
time baseline of high-precision transiting planet surveys 
are the most likely astrophysical false positives. Hence, 
close binary occurrence rates f orm a crucia l input into 
calculations of planet statistics dFressin et al.ll2013 ). 

However, M+M short-period binaries have been enig¬ 
matic because they are relatively difficult to identify vi¬ 
sually due to faintness and spatial resolution limits. Fur¬ 
thermore, predictions for the abundance of such systems 
are inaccessible by state-of-the-art hydrodynamic simu¬ 
lations of star formation, also du e to resolu tion limits 
subjected to computation budget dBatd 120091. [2012 1. To 
date, all empirical constraints o n the statistics of these 
systems der i ve from RV analyses (Fischer fe Marcvlll992l: 
Udrv et al.l l2000t iDelfosse et al.32004 iShkolnik et al l 


Udrv et au ^uuu ; meiiosse et an A 

20101 : IClark. Blake. fc Knappll2012ir 


FM92 find precisely 1 M-dwarf companion within 0.4 
AU among the 62 M-dwarfs subjected to their search. 
Correcting for completeness, they conclude that, in this 
range of orbital distance, the occurence rate of M-dwarf 
binaries is 1.8% ± 1.8%. 

3 Hence the secondary must either be an M-dwarf or substellar. 


Recent work by IClark. Blake, fe Knapd (I2012D (here¬ 
after CBK12) provides an alternative measurement. 
They examine spectra from Sloan Digital Sky Survey 
(SDSS) M-dwarfs to look for RV variations that reveal 
the presence of close stellar companions, thereby deriving 
an estimate on the binarity for orbits with a < 0.4 AU. 
They report the detection of 22 binary candidates among 
the 1452 M-dwarf targets examined. Correcting for de¬ 
tectability, CBK12 infers a close binary fraction of 3-4% 
among the M-dwarfs. 

The NASA Kepler Mission has been highly success - 
ful at advancing exoplanet science (iBorucki et al.|[20Toi ). 
Throughout its primary mission baseline of 4 years, it 
collected continuous and precise broadband photometry 
from some 160,000 stars in a 100-square-degree-patch of 
the sky near Cygnus (field centre: RA = 19h22m40s, 
Dec = +44°30'00"). The resultant light curve database 
has enabled the discovery of thousands of transiting 
planets. Furthermore, the variability in these light 
curves encode information on the activity, rotation, pul- 
sati on, and other characteristics of the stars themselves 
(e.g. Bas tien et aT1l2013HHuber et H1l20l4l Mosser et al.1 


l2014iljohnson et al.ll20l4 lAngus et al.ll2015l ). extending 

their utility beyond those for whom planetary systems 
are fortuitously inclined. Thus, while Kepler may be 
conceived chiefly as an exoplanet mission, it has hailed a 
veritable golden age of stellar astrophysics. 

Eclipsing binary (EB) systems are a class of objects 
observed by Kepler that lie at the intersection of tran¬ 
sits and stellar astrophysics. Furthermore, the occur¬ 
rence of EBs in a sample of targets directly reflects the 
underlying binarity rate, discounted by the geometric 
eclipse probability—a function chiefly of semimajor axis, 
stellar radii, and orbital eccentricity. The detection of 
EBs is a simplified version of exoplanet hunting—a stel¬ 
lar eclipse, even grazing, is generally a much larger sig¬ 
nal than that of a transiting planet, and thus far less 
likely to be missed. In addition to its impressive list 
of planet detectio ns, Kepler has uncovered thousands of 
eclipsing binaries (iCoughl in et al.ll 2011l:lPrsa et al.lf201 lb 

iSlawson et, al.l 120111 : iMatiievic et all 120 12lh Given the 

importance of double-lined spectroscopic EBs as the only 
model-free way of simultaneously measuring the precise 
radius and mass of stars, this is a valuable database 
in and of itself for testing stellar models, especially in 
the regime of low-mass stars wher e theory and observa¬ 
tion are currentl y discrepant (e.g. iLopez-Morales 1120081 : 


iTorres. Andersen, fc GimenezlhoRlT 

Owing to the mission mandate of Kepler , which 
aims to quantify the prevalence of Earth-like planets, 
the methodology for computing exoplanet occ urrence 
statis ti cs has becom e highl y developed (e.g. lYoudjn 


2011 


iHoward ct al.1 12 0121: iDressing fc Charb onnca u 


iPetigura Howard, fc Marcv 


20131 _ iHetig ura. 1__ _ 

Fressin et al.l 120131: iMorton fc Swift] 


Foreman-M ackcv. Hogg , fc Morto n 2014). Since stellar 


systems are akin to planet systems, the same framework 
can be transplanted to calculate stellar companion 
statistics, with several simplifications. The chief source 
of incompleteness ought to be due to non-transiting 
geometries. Since we can quantify the transit probabil¬ 
ity of a binary system given its scaled semimajor axis 
{a/(R\ + 1?2)), which is measured by Kepler , each EB 
detection can be assigned a statistical weight to reveal 
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the underlying population that is not in an eclipsing 
configuration. 

In this work, we explore what Kepler can reveal about 
close (but detached) M+M-dwarf binarity statistics. We 
begin by defining our sample and describing the raw Ke¬ 
pler data and how we conditioned them to our purpose 
(Section [2 ]). T he scheme for culling the EBs is described 
in Section l3Jl The EB detections among the sample are 
listed in Section lT2l and com parisons are made with pre¬ 
vious studies in Section 13.31 Then we elaborate on the 
statistic al m ethod alluded to in the previous paragraph 
(Section [3.41) . The overall NSPS - Number of Stars per 
Star - and its orbital period distribution are presented in 
Section U followed by a discussion of the caveats, biases, 
and comparisons to previous studies. Section[5]concludes 
with a summary of the findings. 

2. SAMPLE & DATA 
2.1. Sample 

Our sample comprises the set of M-dwarf 
targets observed by Kepler as identified by 
iDressing fc Charbonneaul (I2013L henceforth DC 13). 
DC13 also derive updated stellar parameters, such as 
masses and radii, estimated by comparing colours of 
model stars from Dartmouth stellar evolutionary tracks 
with that of the observed stars. For cool stars, these 
parameters should represent a vast impr ovem en t upo n 
the original Kepler Input Catalog (KIC, iBrownl [20111) . 
DC13 provides a total sample of 3905 early to mid 
M-dwarfs with T e g ^ 4000K, each observed for at least 
one Kepler quarter. Distributions of stellar masses are 
shown in Figure [T] exhibiting a heavy but expected 
Malmquist bias toward earlier spectral types, since the 
luminosity-mass relation is very steep (L ~ M 5 ) in this 
region of the H-R diagram. 

Note that the Kepler sample is neither volume-limited 
nor, strictly speaking, magnitude-limited (see Figure [2] 
for the Sloan-r magnitude distribution of the M-dwarf 
sample). As a result, the target selection is poorly de¬ 
fined since M-dwar fs were not the prim ary targets of 
the Kepler Mission dBatalha et al.ll2bl0af) . Nevertheless, 
we argue that this sample is not unrepresentative of M- 
dwarfs in the Solar Neighborhood, particularly for early- 
M spectral types. 

The stars in the Kepler M-dwarf sample did not have 
multi-epoch spectroscopic data at the time of selection, 
and as a result spectroscopic binaries were not intention¬ 
ally excluded. Further, while distances to some stars in 
the Kepler sample extend out to beyond a kpc, the M- 
dwarfs predominantly reside within ~ 200 pc. Thus, we 
expect the occurrence of stellar companions measured by 
the rate of EBs in Kepler to be representative of the bi¬ 
narity for M-dwarfs in the Solar Neighborhood (i.e. for 
relatively nearby stars within < 20 pc of the Sun). 

2.2. Data and Its Processing 

The majority of Kepler’s targets were observed in long- 
cadence (LC) mode, where each data point forming the 
light curve is an integration over « 29.4 minutes. We 
downloaded the Pre-search Data Condition ing S imple 
Aperture Photom etry (PDCSAP) (|Stumpe et al.l [20121 : 
iSmith et al.l l2012h for each target in the DC13 sample. 
This is a publicly available data product supplied by 
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Fig. 1.— The stellar mass distribution of DC13 targets. The 
black solid line depicts all the stars compiled by the DC 13 cata¬ 
logue, of which there are 3905. This means they were observed for 
at least 1 quarter by Kepler. Those stars with at least 3 quarters 
of LC observation, which is a requirement for stars used in our 
statistical study, are shown in red dashes. There are 2975 of them. 
Not surprisingly, the overall sample comprises predominantly early 
M-dwarfs. 
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Fig. 2. — The Sloan-r magnitude distribution of DC13 targets. 
Colour scheme follows that in Figure [T| The sample is neither 
volume- nor magnitude-limited. 


All DC13 stars 

at least 3 quarters of LC 



the Kepler Team and made accessible via the Mikulski 
Archive for Space Telescopes (MAST0). In addition to 
the standard bias and dark subtractions, field flattening, 
and cosmic ray treatment, PDCSAP strives to eliminate 
the systematic and instrumental effects with “cotrend¬ 
ing vectors,” while preserving variations of astrophysical 
origin, such as rotational modulation of spots. 

To prepare the light curves for eclipse searches we per¬ 
form a few additional tasks: we quarter-stitch following 
median-normalization and “flatten” the data over an ap¬ 
propriate timescale using piecewise polynomials. Over 
a chosen window of a set time interval, we perform a 

4 https://archive.stsci.edu/kepler/ 







































4 


low-order polynomial fit to the light curve within this 
window, less any points that are > 5cr away from this 
functional fit, a being the standard deviation of the en¬ 
tire normalized light curve. This is iterated 5 times to 
reach a final polynomial fit for that window. Then the 
window slides over to the next segment and the process 
is repeated, though continuity on the shared boundary of 
the windows is not demanded. The resultant piecewise 
polynomial is subtracted from the original light curve. 
For the chosen polynomial order of 2 and window size 
of 1 day, we filter out smoothly varying features span¬ 
ning longer than this timescale, representing a high-pass 
filtering process. This step rids typical periods of ro¬ 
tational modulation (fast-rotators excepted), while pre¬ 
serving any reasonable detached eclipse event lasting at 
most a few hours for the orbital period range of interest. 
For an early M+M EB on a 90-day circular orbit, the 
transit duration is ~8 hours, and less for shorter orbital 
periods. 

3. METHODS 

Of the 3905 M-dwarfs present in the DC13 sample, 
2975 of them were monitored for at least 3 quarters by 
Kepler , making them uniformly sensitive to sub-90-day 
light curve periodicities. With this sample of 2975 pro¬ 
cessed M-dwarf light curves, we can identify the eclipsing 
binaries among them. The eclipse depth of a stellar bi¬ 
nary is typically much larger than that of a planet transit. 
Of course, in calculating the expected eclipse depth for a 
binary system, one must not only be concerned with the 
radius ratio between the bodies, as one would in a plan¬ 
etary eclipse problem, but also the self-luminous nature 
of the eclipsing body. Hence, the radius and temperature 
of both components matter. 

We generate distributions of the expected primary 
eclipse depths for a realistic population of M+M EBs 
with DC13 M-dwarfs as primaries, accounting for grazing 
scenarios. To each primary we assign secondary masses 
according to an Opik’s mass-ratio power law f(q) ~ q 1 , 
where q = M s /M p , the ratio of the secondary to pri¬ 
mary masses. Then, stellar masses are converted into 
radii and tem perature using observ ationally calibrated 
relations from IBovaiian et al.l (120121) . Orbits are drawn 
from orbital element distributions as described in Sec¬ 
tion 13.41 For M-dwarf binaries below 5 AU, -y may be as 
high as 2.7, but is highly uncertain (jDuchene fc Krausl 
f2013h . We adopt a more modest mass distribution power 
of 1 since it would result in a conservative eclipse depth 
distribution and note that the outcomes presented be¬ 
low are robust to any reasonable variation in 7 . 

The distribution of eclipse depths is shown in Figure 
m Based on our simulations, the eclipse depth can be 
up to 50% for equal-mass constituents, and in 95% of the 
cases the depths are greater than 1.5%, even in grazing 
geometries. Depths of this magnitude are easily identi¬ 
fied in the light curves of even the faintest stars in our 

5 A higher 7 corresponds to greater occurrence of higher-mass 
secondary companions, hence would give rise to larger eclipse 
depths overall 

6 In all the computations, limb darkening and effects 
such as beaming and ellipsoidal variations have been ignored 
(IFaigler &; Mazehj|201U) . Also, we do not consider substellar com¬ 
panions since they are intrinsically rar e aroun d stars generally, and 
in our sample specifically. See Section 14.2.31 


sampliQ is ~ 10. In contrast, the transiting depths of 
most planets discovered to date around cool dwarfs fall 
well below this threshold. 



grazing eclipse depths 

Fig. 3.— A distribution of primary eclipse depths (i.e. fractional 
light blocked in eclipse) for a simulated population of M-dwarf EBs. 
Vertical axis is given in arbitrary units. The red dashed vertical 
line at 0.015 is the cut in depth we make to generate our first 
round of candidates . The magenta da sh-dot histogram shows cool- 
KOI transit depths IlSwift et al.l2015f ). approximated by (fip/f?*) 2 - 
Less than 5 % of the EB cases lie below the selection criterion, con¬ 
trasted with nearly all of the transiting planets around M-dwarfs. 

The eclipse duration is also an important parameter 
that determines the transit window and our ability to re¬ 
solve the shape of the eclipse event given the cadence and 
time baseline of the Kepler observations. The fraction of 
a light curve in eclipse, or the duty cycle, for M+M EBs is 
a function of the total radius, total mass, orbital period, 
and orbital eccentricity of the system. Figure [4] shows 
the distribution of duty cycles for the same population 
of M+M EBs as in Figure [H including grazing geome¬ 
tries. The orbital periods distribution is assumed to be 
a power law of the form /[log(P 0 rb)] ~ [log^orb)]^- For 
the particular simulation shown in Figure [H p = 2. The 
period is restricted to below 90 days, or what is feasible 
for an unambiguous detection and periodicity determina¬ 
tion with at least 3 quarters of Kepler data. All but 0.3% 
of the duty cycles are greater than 0.04%. This outcome 
varies little over a wide range of reasonable choices for p. 

Given the eclipse depths and durations predicted by 
our simulations we can devise the relevant criteria for re¬ 
liably identifying EB candidates with an automated pro¬ 
cedure. After the automated process identifies potential 
EBs, we then investigate the light curves by eye to iden¬ 
tify the final sample, as we describe in further detail in 
the following section 

3.1. Identifying the M+M-EBs 

We confine our search to all stars with at least 3 quar¬ 
ters of Kepler observations (of which there are 2975 ob¬ 
jects) such that we are uniformly sensitive to EBs with 
periodicities of < 90 days. We use an automated al¬ 
gorithm to exclude the obvious non-EBs from the sam¬ 
ple, while retaining those with suggestive features of EB 
candidates. We then run a semi-manual visual vetting 

7 Typical S/N for a 1.5% single-point transit for M-dwarfs with 
Kepler magnitude (K p ) > 16 
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Eclipse Duty Cycle (T dul /P orb ) 


Fig. 4.— A distribution of theoretical total primary eclipse dura¬ 
tions as fraction of orbital period for a general simulated M-dwarf 
EB population in the same manner as for Figure [3] Vertical axes 
are arbitrary units. The red dashed line marks the position of our 
first-round duty-cycle cut on the light curves. It lies below 99.7% 
of the eclipsing cases, missing only the very grazing ones. Note 
that duration of ingress/egress and finite cadence would reduce 
the number of points expected to be captured at the theoretical 
maximum primary eclipse depth. 

routine on the relatively small number of targets that 
survive the initial vetting. Each step is described below. 

3.1.1. The Depth Cut 

Our automated search algorithm first checks to see if 
each light curve has downward flux variations that meet 
the depth criterion informed by our simulations. Figures 
[3] and [4] demonstrate that most M+M EBs have primary 
eclipse depths above 1.5% and duty cycles above 0.04%. 
For a Kepler quarter of ~90 days at ~30 min cadence 
containing ~4300 data points (neglecting intra-quarter 
data gaps), this amounts to requiring 2 or more of the 
points in the conditioned and normalized light curves to 
be below 0.985. Using this criterion, we end up with 150 
candidates, which we subject to visual inspection. 

3.1.2. Semi-Manual Inspection in Frequency Space 

For the remaining candidates, we perform a series of 
Fourier transforms on the un-flattened, i.e. PDCSAP, 
timeseries in a supervised hierarchical manner to decide 
whether a periodicity exists and, if so, whether it is as- 
cribable to stars eclipsing one another. This is done 
through a custom user interface written in IDL. The de¬ 
cision to use the PDCSAP rather than the flattened light 
curves from the previous step is to allow the human to 
assess all astrophysical features holistically and to better 
identify the most likely cause of the feature. 

Interpolating the data to even spacing in time (i.e. 
at 29.4 min intervals), we then perform a Fast Fourier 
Transform (FFT), resulting in a power spectrum. When 
interpreted effectively, the FFT procedure can reveal 
much informatio n about the causes of periodicity in a 
time series (e.g. iRannanort et al.ll2014l h A light curve 
with no periodic variation will contain no significant 
peaks in the FFT. A quasi-sinusoidal signal with a wide 
duty cycle, as that associated with rotational modula¬ 
tion or BEER effects (BEaming, Ellipsoidal, Reflection - 
brightness mod ulations associated wit h binaries but not 
their eclipses) dFaigler fc Mazehl 1201 111 , generally shows 


up as a prominent and narrow fundamental peak fol¬ 
lowed by a rapid decay of harmonic peaks. Repeated 
narrow features, such as the signature of an EB, results 
in many repeated peaks (i.e. fundamental + harmonics), 
perhaps alternating in height if the secondary eclipses are 
sufficiently prominent, with the overall height structure 
modulated by an envelope that is only slowly declining, 
(see bottom panel of Figure [5]). The central peak fre¬ 
quency usually does not exactly equal, but is very nearly 
the true periodicity or an integer fraction multiple of it. 
The contrast with smoothly-varying periodicities is the 
inherent nature of frequency transforms: narrow features 
in the time domain translates into broad features in the 
frequency domain, and vice versa. 

When significant peaks are detected in the Fourier 
transform, they are immediately apparent to the eye. 
It turns out that most of the candidates exhibit some 
apparent peak structure in their FFT spectrum, though 
most of these appear to be due to the rotational modu¬ 
lation of large-amplitude fast-rotators as opposed to pe¬ 
riodic eclipses. Such objects have escaped our vetting 
process unscathed because their deep rotational modula¬ 
tion is also too rapidly varying to be fitted out without 
obscuring the integrity of possible eclipse signals. 

As a final step, we phase-fold the light curves of po¬ 
tential EBs to search for the unmistakable V-shape at a 
very precise periodicity. In the phase-folded light curve, 
a non-EB might show large-scale variations akin to rota¬ 
tional modulation, as well as some phase-drifting, while 
the EB would already exhibits an obvious dip structure. 
Further period refinement results in a light curve akin to 
that displayed in the top panel of [5] 

When refining the true EB periodicity for a 
promising target, we also perform a Box-Least 
Squares (BLS) search around the suspected period 
dKovacs. Zucker. fc Mazehll2002l h This results in a more 
precise periodicity, though still not perfect, for the BLS 
was designed to find planet-transit-shaped events with¬ 
out obvious secondary eclipses. The final outcomes are 
EB orbits determined to a suitable precision for our work 
(i.e. typically ~ 10“ 3 days), and that we have in princi¬ 
ple found all the unambiguous EBs within some period 
range to which the Kepler observation baseline ought to 
be sensitive. 

Running this process on our sample of M-dwarf light 
curves results in 17 EB candidates with periods from 0.3 
to 76 days. Since we are investigating the occurrence 
of unambiguously detached systems whose periods are 
greater than 1 day, we discard the 3 very short-period (< 
1 day) targets, namely KID4936334, KID9077796, and 
KID12004834, which exhibit various degrees of contact. 
This is evident in the duration of their transits, which 
much exceed the upper limit of expected M+M detached 
EB transit durations given their orbital period. Their 
phase-folded light curves are shown in Figure [9] This 
criterion whittles our candidate list down to 14. 

3.2. Further Diagnostics and the Final List 

The EB candidates are examined further via a couple 
of simple diagnostics. We inspect the eclipse duration as 
well as centroid-flux correlations. The former test may 
indicate whether the eclipse could be associated with a gi¬ 
ant. The latter validation technique is standard in iden¬ 
tifying blends and false positives for planet detections 
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Fig. 5. — Top panel: A sample flattened and phased light curve 
of the EB KID3830820. Middle panel: its Fast-Fourier Transform, 
plotted in terms of period = 1/frequency. The range of periods 
plotted is 0.1 to 200 days. The red dot, pointing to the longest- 
period prominent peak, is the true period of the system upon which 
the light curve is folded. Bottom Panel: the FFT plotted in terms 
of frequency. Note the slowly declining alternating peak structure, 
a signature of EBs. 

(e.g. iBatalha et al.ll2010bHTorres et all 1201 111 . 

The eclipse duration check involves comparing the ap¬ 
parent duration of the eclipse with that expected for 
the fiducial stellar parameters. If the actual duration 
much exceeds this limit, then the eclipse may be asso¬ 
ciated with a giant star or a very eccentric orbit, or is 
inconsistent with the detached EB interpretation. While 
the 3 sub-day candidates do exhibit eclipses of unusu¬ 
ally long d urations and shapes suggestive of contact (Sec¬ 
tion [37TT2]) , the remaining 14 EB candidates show eclipse 
durations consistent with their orbital periods and esti¬ 
mated stellar radii. 

We also compute the positions of the brightness cen¬ 
troids in the pixel-level ‘postage stamp’ data using a 2D 
Gaussian fit. Correlating the centroid position with flux 
could show photocentre shifts throughout the light curve. 
A significant shift during eclipse would suggest blending 
and light contamination, raising the possibility that the 
source of the light loss does not belong to the foreground 
target centred in the aperture. 

One candidate shows dramatic centroid shifts during 
eclipse. KID3830820 is blended with two other stars 
(~ 3" and ~ 20"), both of which fall partially into the 
main aperture. The stars may be physically associated. 
Pixel-level light curves suggest that the eclipses are in¬ 
deed associated with the main M-dwarf targeted. 

Two more candidates out of the original 17, 
KID5794240 and KID11548140, have survived all our 
search criteria but follow-up data have revealed that 


their eclipsing companions are non-stellar. KID 5794220, 
or KQI-254.01, is a n M-dwarf hosting a giant planet 
(| Johnson et alllMll . KID11548140, or KOI-256.01, is 
a post-com mon envelope bi n ary w hose companion is a 
white dwarf dMuirhead et al.ll2013l l. 

Our final 12 M+M-dwarf EB candidates are listed in 
Table [TJ None of these stars started as Guest Observer 
targets. Their phase-folded light curves are displayed in 
Figure [7] and Figure |U Note that all the light curves 
show V-shaped eclipses and/or prominent secondaries, 
characteristic of EBs. 

The longest orbital period is 76 days, or 0.21 years (log 
-Porb = -0.68), and many are below 3 days, or 0.008 years 
(log P or b = -2.1). The fact that we have few detections 
at long periods and many in close orbits is consistent 
with expectation despite ignorance of the actual binarity 
distribution, for the average geometric transit probabil¬ 
ity declines rapidly with orbital period, and the growing 
likelihood of non-zero eccentricities at larger orbital sep¬ 
arations make it less likely to observe both the primary 
and secondary. Hence we would expect fewer wide-orbit 
EBs to be detected. This expectation is an observational 
bias that relates to the fact that longer period EB detec¬ 
tions have the greatest influence in our conclusions about 
overall binarity rates, for they carry the greatest statist- 
cal weight, being the least likely to be found. We discuss 
the statistical framework from which we draw inferences 
in Section [TU 

3.3. Comparison to existing Kepler EB catalogues 

EBs in K epler have been tracked and studied by previ¬ 
ous groups. iCoughlin et al.l (1201111 catalogued and mod¬ 
elled the light curves of all detached EBs with primary T 
< 5500K and orbital periods under 32 days in campaign 
1 of Kepler. The Kepler Eclipsing Binary Working Group 
have produced the Villanova catalogue, a running compi¬ 
lation of all detected EBs in the Kepler datase t (for meth¬ 
ods and a description of t he catalogue, see iPrsa et al.1 
120111 : iSlawson et al.l 120111 : iMatiievib et al.1I2012H . The 
catalogue contains classification of each candidate, pro¬ 
viding a benchmark against which we can compare the 
outcome of our search procedure, and to cross-check their 
catalogue completeness for M+M EBs. Within the DC13 
sample, there are 19 overla ps with the Villanova Cata- 
logue, which includes all the ICoughlin et al.l (l201lh EBs 
in the same mass range. All seventeen of our original 
EB candidates are among this overlap. Of which three 
(KID4936334, KID9077796, and KID12004834) satisfy 
our two initial cuts but possess orbital periods within 0.4 
days and unusually large duty cycles, making their inter¬ 
pretation as normal detached binaries unreliable. They 
have therefore been excluded and we have limited the 
scope of our study to periods above 1 day. These very 
short-period binaries are shown in figure [9] Yet another 
two ‘EBs’, KID5820218 and KID9892651, are unconvinc¬ 
ing upon visual inspection. Figure [TO] shows their re¬ 
spective light curves folded upon their documented peri¬ 
ods. The eclipses, if actually present and due to stellar 
sources, are very low signal-to-noise and would represent 
either a marginally grazing geometry or a nearly substel- 
lar eclipsing companion, or both. 

We have encountered a few other unambiguously 
detached M-dwarf EBs not documented in the Vil¬ 
lanova Catalogue, namely KID5769943 (P or b = 1-01 d), 
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TABLE 1 

Final List of M-dwarf M-dwarf Detached Eclipsing Binaries in the Kepler Field 


KID 

-^"orb 

(days) 

s* . 

prim 

Tdur /-^orb 

K p 

(mag) 

# Quarters 
Observed 

M* 

(M G ) 

R * 

(Rq) 

CDPP-3 ** 
(ppm) 

Eccentric? *** 

Notes 

2442084 

49.789 

0.251 

0.00603 

15.59 

17 

0.548 

0.541 

362.24 

Yes 

1 

3830820 

15.583 

0.063 

0.0104 

15.37 

17 

0.563 

0.555 

314.47 

Yes 

2 

5871918 

12.643 

0.256 

0.0173 

15.70 

17 

0.573 

0.562 

317.21 

Yes 

3 

6023859 

27.010 

0.079 

0.00156 

15.48 

8 

0.559 

0.537 

368.92 

Yes 

4 

6620003 

3.429 

0.033 

0.0271 

15.686 

17 

0.542 

0.534 

287.35 

No 

5 

7605600 

3.326 

0.193 

0.0296 

14.888 

8 

0.359 

0.350 

376.70 

No 

6 

7671594 

1.410 

0.106 

0.0355 

15.815 

17 

0.245 

0.250 

3305.8 

No 

7 

9772531 

31.202 

0.040 

0.00319 

15.798 

17 

0.462 

0.450 

281.53 

Yes 

8 

10979716 

10.684 

0.091 

0.0135 

15.774 

17 

0.553 

0.524 

268.35 

Yes 

9 

11546211 

2.194 

0.041 

0.0301 

15.155 

13 (inch Ql, Q17) 

0.339 

0.334 

488.50 

No 

10 

11853130 

76.87 

0.017 

0.00135 

15.949 

13 (inch Ql, Q17) 

0.407 

0.380 

393.58 

Yes 

11 

12599700 

1.018 

0.051 

0.0518 

15.780 

17 

0.493 

0.480 

1106.2 

No 

12 


Notes: 

*: Sprim is the approximate depth of the primary eclipse. 

**: CDPP-3 is the ‘Com bined Differential Photom etric Precision’ over 3 hours, a metric for the noise in the light curve of the star over a 
typical transit timescale d Christiansen et al.H20l3) . Here we displayed the mean CDPP-3 over all observed quarters. 

***: Eccentricities are judged from the phase difference separating the primary and secondary eclipse. A phase separation of 0.5 is 
consistent with a circular orbit, provided the argument of periapse, uj, is not exactly 7r/2 or 37 t/2 . For other values of phase separation we 
assume the orbit to be eccentric. This includes systems for which only the primary is visible. 

1: unambiguous and classic primary and secondary eclipses. Phase separation (0.181) suggests considerable eccentricity. 

2: deeper primary accompanied by very shallow secondary with phase separation of 0.565 suggesting non-circular orbit. This star is 
blended with two stars (visual separation ~ 3" and ~ 20" whose lights contaminate that of the chief target to differing extents from 
quarter to quarter. As a result, this target exhibits dramatic centroid shift in eclipse. Pixel-level data suggest that the eclipses are indeed 
associated with the main M-dwarf target, though the depths ought to be treated with caution, since blending causes considerable variation 
in absolute and relative brightness from eclipse to eclipse. 

3: clearly visible rotational modulation signature on top of the eclipse signals. The eclipses are not aligned with the rotation phase. Orbit 
is eccentric. 

4: the eclipse is suspiciously sharp, i.e. short-duration, perhaps betraying the eclipse being rather marginally grazing. The true 
orbital period could be twice that recorded here, for visual inspection reveals the transit depths are subtlely different if folded upon 
P orb = 2 x 27.010 = 54.020 days. In fact, RV data suggest this may be the case (Jonathan Swift, private communications). 

5: distinct primary and secondary eclipses. Clear rotational modulation signature, appearing to be synchronized with the orbit, which has 
also been circularized. 

6 : notable rotational modulation in amplitude, with one periodicity being 3.9 days. Occasional but sizable flares between periods of 
quiescence. Orbit is consistent with circularity. 

7: strong evidence of rotation-orbital synchronization apparent in minimal scatter of the underlying rotationally-modulated phase curve 
when the light curve is folded upon the orbital period. 

8 : very sharp eclipses on a very quiet star. Only primary is visible. Rounded bottom suggests the eclipse geometry is not grazing. 
Eclipse depth indicates the eclipsing body is ~ 0.09 -Rq, which does not rule out brown dwarf or Jupiter-sized planet as the companion. If 
interpretation as stellar binary were robust, then the orbit must be fairly eccentric for the secondary to be out-of-view. 

9: clearly visible rotational modulation. Orbit is not fully circularized. Effect of eclipse timing variation (ETV) is appare nt, as it smears 
out the period-folded eclipses. A third, farther body in the system inferred from ETV has been reported (jBorkovits et al1l2Q15IL 
10: flagged as ‘false positive’ in MAST and catalogued as an EB in Villanova. Large-amplitude rotational modulations with periodicities 
2.06 and 2.16 days, strongly suggesting rotational synchronization with the orbit. 

11: the longest-period EB in this sample. Only primary is visible, and it is shallow and narrow. If stellar in nature, these characteristics 
suggest the eclipse geometry is grazing with eccentr i city. It is also not inconsistent with interpretation as a giant planet (e.g. in a 
most recent planet candidacy catalogue, I Swift et all (120150 includes this target as a cool-KOI, KOI-3263.01, with listed radius being 
0.057-R®, albeit with a hig h false positive probability). Lucky Imaging with AstraLux at Calar Alto Observatory (PI Lillo-Box, 
(|Lillo-Box et al.1120121 . [20140 j and recent adaptive optics with NIRC2 at Keck Observatory (PI Ciardi) reveal a faint visual companion 
at '•'•'0.8”. Followup observations are continually updated through the Kepler Community Follow-up Observation Program, CFOP 
(https://cfop. ipac.caltech.edu). 

12: a restless and regularly flaring system. Notable rotational modulation synchronized with eclipses. 


KID7938883 (P orb = 0.490 d), KID8949316 (P orb = 
0.604 d), and KID12009213 (P orb = 19.4 d). However, 
they are all part of the Guest Observer (GO) program 
GO2OOO10, whose targets were observed for one quarter 
only (between Q6 and Q9), hence do not fulfill the cri¬ 
teria of having an observing baseline of at least three 
quarters. 


8 Interestingly, G020001 contains 1196 targets and was intended 
to identify new short-period M-dwarf eclipsing binaries for stellar 
characterization. 


3.4. Statistical Analysis and the NSPS 

Since Kepler’s launch, a number of statistical studies 
have been published to infer the occurrence and dis¬ 
tribution of planets around their host stars, studying 
dependencies on planet radius, mass, orbital period, 
eccentri city, planetary mult i plicity, and orbita l align¬ 
ment j[Howa^^nd]_[2(^2L_[SHisinf^E-Cha^onneau 


2013 : iPetigura. Howard, fc Marcvl 120131 Fressin et al 


20131 iForem an-Mackey. Hogg, fe Mo rton _ 12014 

Ballard fe Johnsonl 120151 : [Morton fe Swiftl 12014 here¬ 
after MS14). Most of these studies provide a measure 


































the number of pl anets per star , also known as NPPS 
(see, e.g., MS14). lYoudinl (|2011l f provides a generalized 
statistical framework for measuring NPPS using a max¬ 
imum likelihood analysis assuming Poisson statistics. 
The non-parametric version proceeds as follows to give a 
point-estimate of companion occurrence integrated over 
the entire phase space of system properties searched. 
Referring to the definition of NPPS in their Equation 9: 


fi 


Ni 

Vi N *’ 


(1) 


where TV* is the total number of stars in the survey, l 
labels one of the multiple bins into which detections are 
divided (e.g. a range of orbital periods), Ni is the num¬ 
ber of planet detections in bin l , and rji is the detection 
efficiency in that bin. In effect, 1/r/i serves as a ‘statisti¬ 
cal weight’ on the fraction of planets detected within the 
survey star sample, which corrects the actual detections 
to the true number that exists inside the bin. Summing 
fi over all bins give the overall planet occurrence: 


NPPS = ^/ ; . (2) 

i 


The problem of stellar multiplicity far predates that 
of planet multiplicity (see Section [TJ, but the techniques 
used to measure planet occurrence from a transit survey 
are directly analogous to our goal of measuring binarity 
from stellar eclipses. The above methods constitute the 
starting point of our analysis which, in many respects, 
is a simpler problem than exoplanet statistics; detection 
of an object as large as a star ‘transiting’ a star is not 
hampered by low signal-to-noise, unlike detecting planet 
transits. Also, since we are concerned only with M+M 
dwarf binaries, several simplifications can be made. As 
has been demonstrated in Section 13.1.11 all but the most 
grazing of eclipsing geometries ought to produce signals 
that far exceed the amplitude of, and have very different 
shapes than typical stellar photometric variability. Thus, 
false positive rates should be low for our study and the 
discovery efficiency is chiefly limited by the geometric 
eclipse probability. 

In general, the eclipse probability for binary stars is: 


Ptr(Rtot,a, e, w) = Max 


Rtot 

a 


lie sin 
1 — e 2 A 


( 3 ) 


where R t ot = R* l + R* 2 , e is the orbital eccentricity, u> 
is the argument of periastron, and a is the semi-major 
axis, related to the orbital period by Kepler’s 3rd Law. 
The ‘A’ in the A refers to the primary, whereas the ’ 
applies to the secondary eclipse. Marginalizing over all 
w € [0 , 27t] results in 


Ptr(R*) M k . Porb) — Rtot 


r 4 tt 2 1 1/3 / 1 \ 
PLGM tot \ U-(e 2 )A 


( 4 ) 

where M tot = + and (e 2 ) is the average squared 

eccentricity of the population of binaries. For main- 
sequence stars, R* « M*, in Solar units, and as a result 
R tot ~ M to t- Eliminating constants, we may write the 


following approximate dependency: 


Pti(Mt 0 t, Rorb) 




Porb 


2/3 



( 5 ) 


More massive systems with shorter orbital periods are 
more likely to transit. A population with a higher aver¬ 
age eccentricity is more likely to eclipse. The distinction 
with exoplanet transits is that both components of the 
binary system contribute. Pt r is taken as the detection 
efficiency r)i, where bin l is identified by the system or¬ 
bital period. 

When dealing with a large sample of stars, the prob¬ 
ability of transit around every star must be accounted 
for. In practice, this condition modifies 77 for each detec¬ 
tion. As MS14 puts it, 77 should be a factor that obeys 
the following statement: ‘if a very large number of stars 
identical to star i were distributed randomly around all 
the stars in the survey, only a fraction rji could have 
been detected.’ We treat this issue by essentially aver¬ 
aging over the theoretical 77 around each primary target 
star in the survey sample. 

In the small-bin-limit, we can compute the Number of 
Stars Per Star , NSPS(M*,P orb ) using the correspond¬ 
ing 77(M*,P or b) to estimate the number of undetected 
systems based on each detection. Binning the results 
into ranges in parameter space of interest, be it in M* 
or P or b, we may obtain the distribution function in these 
chosen parameters. This would be the case if the system 
configurations are perfectly known. 

However, an eclipse light curve alone, unsupplemented 
by spectroscopic information, encodes a limited amount 
of information. While orbital periods can be measured 
very accurately, it is difficult to estimate the primary 
and secondary masses of each given system. So while 
vast simplifications occur for EB detections, the non- 
negligible size and self-luminous nature of both compo¬ 
nents involved introduces some additional complications. 
To compute 77 (M+, R*, P or b) in a meaningful way for each 
detection, we use a Monte Carlo approach. Each target 
in our sample, be it single or multiple, is attached to a 
mass and a radius value derived from DC13’s isochrone 
fits. The primary mass M*i and radius R*i are equated 
with these values. This is a reasonable assumption be¬ 
cause, for equal-mass systems, the photometry of the sys¬ 
tem is not too dissimlar to each individual star and, for 
low mass-ratios, flux from the primary will dominate the 
system flux. 

We draw the the secondary star’s mass, M* 2 , from a 
mass ratio distribution of form f{q) ~ q , where the ra¬ 
dius is calcula ted using the empiri cal mass-radius rela¬ 
tionship from iBovaiian et all 120121 Next, we draw the 
orbital eccentricity (e) from an eccentricity probability 
distribution and the argument of periastron (w) from a 
uniform distribution, and the inclination (?) from a dis¬ 
tribution uniform in cos(«). 

Unfortunately, the true eccentricity distribution for 
close M+M binaries is not well measured. An exploration 
of this issue by an e xtensive RV survey for low-mass spec- 
troscopic binaries (jMavor et al.i 1 2000 : ilJdrv et al ii 2000 : 
show that binaries with components on or beyond the 
main-sequence and periods < 10 days have been circu¬ 
larized, while for periods longer than 10 days the ec¬ 
centricity distribution is a steeply rising function of or- 
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bital period. iRaghavan et al.l (i2010f) ’s study of solar- 
type stars report similar findings, namely orbital circu¬ 
larization for binaries under 12-day periods. Beyond this 
threshold, any e c centri city up to 0.6 is equally likely. 
iHalbwachs et alJ (|2003h observe that, for FGK bina¬ 
ries, those with higher mass ratios preferentially possess 
dampened eccentricities relative to their more unequal 
mass ratio counterparts. These findings are consistent 
with tidal evolution theory. The tidal circularization 
timescale, which governs how likely a binary system is 
found circulariz ed, is a ris ing function of orbital period 
and mass ratio (Hu mm- 

Despite these clues, no functional forms exist for a joint 
probability distribution of eccentricity in terms of orbital 
period and mass ratio. The marginal cases in our own 
EB sample, namely two systems with P or b ~ 10 days and 
Porb ~ 12 days (refer to Table [T]), are both eccentric, as 
evidenced by the timing of secondary eclipse relative to 
the primary eclipse. Therefore, we adopt a piecewise 
eccentricity distribution, where e (P or b < 10 days) = 0 
and e (P or b > 10 days) £ [ 0 . 0 , 0 . 6 ] uniformly. We also 
perturb the upper bound between 0.4 and 0.8 to test the 
robustness of the final statistic to this assumption. 

The process for simulating orbits for each EB detec¬ 
tion is repeated 1000 times around every star surveyed, 
yielding a distribution in detection efficiency 77 at every 
detected EB orbital period, fj(P). The detection is then 
weighted by ry _1 and normalized by the total number of 
survey stars to reach a distribution of fi as defined in 
Equation ffl. From here, we can give our best point es¬ 
timate of NSPS via eqn j2|). In the small-bin limit where 
each detection i constitutes its own bin, NSPS = )P j: /,. 

Since the actual number of detections is quite small, to 
reach a period distribution of the binaries via the typical 
route of constructing a weighted histogram is problem¬ 
atic. In particular, binning in histograms is arbitrary 
and discontinuous, especially if each bin contains only a 
handful of detections. A subtle change in the binning 
scheme could drastically change the shape of the distri¬ 
bution. 

To avoid this issue, we adopt the maximum likelihood 
analysis for estimating a parametric distribution of com¬ 
panion occurr ence following the procedure detailed in 
lYoudinl (EvTilf . So long as the target distribution can 
be reasonably parametrized, this approach does not in¬ 
volve binning, and also yields confidence intervals. We 
begin by writing a general differential distribution: 

df(x) _ . , . . 

—q— =Cg a (x) (6) 

so that df = {df /dx)dx represents the probability of a 
primary M-dwarf star having a binary companion in the 
interval dx. x is some general property of the system 
and a is a shape parameter that defines the functional 
form of the shape function, g. C is an amplitude or 
normalization factor. 

We set x = log 10 P and a to be a power law expo¬ 
nent that parametrizes the underlying period distribu¬ 
tion, such that: 

g a ( X ) = e ax = (J0 (7) 

here Pq is a reference period value, which is set to 1 day. 


The total number of binary companions around TV* 
surveyed stars is then: 


N+ 

N tot = CJ2 

3 


g a {x)dx 


( 8 ) 


as in eqn (15a) of lYoudinl (120111 1. 

We c an defin e a ‘sh ape integral’, F , analogous to eqn 
(17) of lYoudinl (1201 111 , to weigh the shape function over 
all the stars surveyed and the detection efficiency: 


1 f C 

/ Vj(x)g a (x)dx = / fj(x)g a (x)dx 


(9) 


where 


V(x) 


1 

A* 


iV* 

'52, r h( x ) 

3 


( 10 ) 


this is essentially what we mean by computing fj(P) ear¬ 
lier via Monte Carlo simulation. The expected total num¬ 
ber of binary detections, given g a (x), is N ex = N+CF a . 
rjj is the net detection efficiency of a companion around 
star j with system property x. 

The likelihood of the observed data, i.e. number of 
detected EBs, N eb , is a product of their Poisson likeli¬ 
hoods: 


L = 


Web ' 

n 

. 2=1 


exp (—./Vex) 


( 11 ) 


which, given the formalism introduced thus far, can be 
rewritten as: 


L = 


N e b 

C N ° h W_g a {,Xi) 
2=1 


exp(—A ex ) 


( 12 ) 


with the planet parameter differential, dx , eliminated 
since it represents a constant interval. It turns out that 
the logarithm of th e likeli h ood in eqn[T3]can be put in the 
following form (see lYoudinll201ll and references therein), 
the sum being performed over all EB detections: 


In L a = -N eh lnF a 


N e b 

+ y]ln(ff a (a;j)) + N eh 




(13) 


the best-fit normalization, C, is the one that maximizes 
likelihood, i.e. 9 In L/dC = 0 and is precisely: 


C n 


N e b 
N*F a 


(14) 


i.e. the most likely C is that which matches the expected 
number of detection, N ex and that actually detected, 
N eh . 

Our goal is to find the best-fit a and its confidence 
intervals. We do so by searching over a grid of a’s 
and computing the log-likelihoods of the data as given 
in eqn 1131 calculating g a and F a for each. Contours 
in log-likelihood can be compared to a Gaussian dis¬ 
tribution. Around the maximum likelihood, ln(A max ), 
the na confidence interval involves a values that results 


































10 


in ln(L Q ) ^ ln(L max ) — n/2. The NSPS is given by 
NSPS = Nt 0 t/N+, Ntot from eqn ©. The best-fit value 
is always consistent with that derived from eqn (j2|) . 


4. RESULTS AND DISCUSSION 

We have selected M-dwarfs from the Kepler Mission 
with more than 3 quarters of photometry and identified 
the eclipsing binary systems among them. These consti¬ 
tute all the M+M close binaries with favourable orbital 
inclination geometries to result in detection as EBs. A 
statistical framework akin to that used for exoplanetary 
occurrence rates is employed to analyze the data and 
infer the orbital period distribution and overall close bi- 
narity rate of the detached M-dwarfs in the Kepler held. 
We summarize the result of the exercise described thus 
far with Figure[6] Integrating our derived power-law dis¬ 
tribution over the period range spanning 1-90 days re¬ 
sults in NSPS= 0.11 ± 0.02. Perturbing the eccentricity 
distribution by setting the maximum eccentricity to 0.4 
and 0.8 result in a systematic 0.01 upward and downward 
adjustment to the median, respectively. 



Fig. 6. — The inferred power-law orbital period distribution and 
overall occurrence of M+M binaries in the Kepler field from EB 
detections. The best-fit power-law exponent, a, and its la interval, 
is 0.52 ± 0.22. This corresponds to NSPS = 0.106 ± 0.02. The red 
vertical lines delineate the actual EB detections and their relative 
weights due to detection efficiencies. 


4.1. Previous Studies 

The most directly comparable published studies of M- 
dwarf multiplicity sensitive to this tight period range are 
FM92 and CBK12, both of which use RV variation over¬ 
time as an indicator of stellar companionship within a 
few AU. In the case of FM92, the typical dataset for a 
primary target contains 15 spectra over 4 years. In prin¬ 
ciple, the companion sensitivity is well below the substel- 
lar threshold. They conclude that the binarity fraction 
within 0.4 AU is 1.8 ± 1.8%. CBK12’s study is based on 
RV variability measured from stellar spectra in the Sloan 
Digital Sky Survey (SDSS) over 2-30 days. They infer a 
binarity fraction in a comparable orbital distance to be 
3 — 4%. The present study uses photometric time series 
data and the detection of eclipsing binaries to argue a bi¬ 
narity fraction of 11%±2% within a 90-day period. Both 
previous results are in tension with the present study. 


In what follows we examine the particulars in the pre¬ 
vious binarity studies, from sample selection to analysis 
technique, to highlight the merits and limitations of each, 
and to assess the appropriateness of direct comparison. 

4.1.1. FM92 

The FM92 sample consists of stars targeted as part 
of an RV program to search for substellar companion s 
around mid- and late-M-dwarfs (|Marcv fc Benitzlfl989 1. 
As such, the affiliated spectral identifications of their ob¬ 
jects are all M2 or later with V<11.5 mag. To accen¬ 
tuate the dynamical signature of brown dwarfs, which 
is the original objective of this observing program, the 
FM92 RV sample excluded targets that possess a known 
companion within 10" to avoid two stars falling on the 
spectrometer entrance slit. Since targets amenable to RV 
studies must be very bright, the final sample in FM92 
consists of 62 M-dwarfs in a magnitude-limited sample 
all residing within ~ 15pc of the Sun, at declinations 
of —10° to +50°. The authors assess completeness fac¬ 
tors in terms of orbital elements and companion mass, 
and adopt a conservative estimate of this factor to derive 
the total number of companions. Among these M-dwarf 
primaries, only a single companion was detected within 
0.4 AU. Thus, the measured binarity in the semi-major 
axis bin 0.04 - 0.4 AU is drawn from this sole example. 

FM92 presumably use the central limit of the Poisson 
distribution to arrive at their binary rate of 1.8 ± 1.8%. 
Note that, with only a single detection, it would be more 
appropriate to use binomial statistics, which gives rise to 
an asymmetric probability distribution of 2.8jt^'g%, i.e. 
with a larger tail towards higher values. 

4.1.2. CBK12 

CBK12 use M-dwarf spectra collected from SDSS. The 
sample is nominally magnitude-limited (the initial selec¬ 
tion criterion enforces 16 < i < 20.5), and in total con¬ 
tains 39,543 M-dwarfs (though in practice, the sample 
used for companionship search is much smaller, see be¬ 
low). 

To contrast with the FM92 targets, whose spectra col¬ 
lection is more systematic, the SDSS spectra are recorded 
somewhat haphazardly over very short time baselines 
and the sampling is relatively sparse. The intent of SDSS 
spectra is for classification and characterization as op¬ 
posed to monitoring for radial velocity variations. In 
other words, it is a serendipitous RV survey. 

In their analysis, CBK12 used a combined spectrum 
formed by co-adding all exposures for each star. Typi¬ 
cally they are three 15-minute exposures taken sequen¬ 
tially, such that the time baseline spanned lies within an 
hour. In rare circumstances, though, the set of spectra 
used to construct the combined spectra would be sep¬ 
arated by baselines as long as days, lending themselves 
to the possibility of observing radial velocity variations. 
Within the M-dwarf sample, 1452 stars have observa¬ 
tions spanning over 2-30 days. It is within this M-dwarf 
subsample that the authors searched for RV variations 
indicative of companions. 

The SDSS M-dwarfs are large in number and span 
a relatively large volume in space. The stars are sam¬ 
pled in small patches that are u niformly distributed in 
the northern sky (5 > —20°) (iStoughton et al.1 120021 : 
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lAbazaiian et al.l 120091 lYannv et al.l I2009H . However, 
CBK12 is not volume-limited and not really magnitude- 
limited, either. Referring to Figure 1 therein, the ac¬ 
tual magnitude distribution of the sample is unimodally 
peaked at i « 16.8 and tapers to nearly none at * ~ 18. 

Many SDSS stellar spectra have S/N > 30, sensitive 
to radial velocities down to ~ 4km/s, making them 
amenable to identifying bulk kinematic properties of 
galac tic structure, such as stellar streams (lYannv et al.l 
120091 ). However, they are not ideally set up for RV work. 
The binary detections are identified by assuming all de¬ 
tected RV variability above a threshold can be attributed 
to the influence of a binary companion. The details of 
the criterion itself assumes an underlying and station¬ 
ary error distribution of RV variability, which is inferred 
from a control sample of targets whose RV time baselines 
are < 2 hours. Applying this criterion, there are 22 de¬ 
tections throughout the CBK12 sample. After complete¬ 
ness corrections, they conclude a close binary fraction 
of either 3% or 4% within 0.4 AU, depending on their 
assumption of the frequency distribution in semi-major 
axis. The Bayesian posterior la confidence interval is 
~ 1% for each scenario. 

4.2. Particular Points of Comparison 
4.2.1. Sample Comparability 

Are we really all probing the same population of ob¬ 
jects? As mentioned under Section 14.1.11 the primary 
stars in the FM92 RV sample are predominantly of spec¬ 
tral type M2 and later. In contrast to the typical Kepler 
M-dwarf target, which lies ~ 200pc away, the FM92 stars 
reside in the very local solar neighbourhood (<15pc). In 
addition, all potential targets harbouring a known com¬ 
panion within 10" have been excluded. It is now well- 
known that spectroscopic and con tact binaries are li kely 
to have distant triple comp anions dTokovinin et al. 1120061 
iPribulla fc Ruc inski 200(f), whereas FM92 assumes this 
correlation to be negligible. Hence the cut made by 
FM92 could bias their sample in the search for close- 
in companions. The FM92 sample is magnitude-limited 
(V < 11.5), though the targets are ultimately chosen for 
their amenability to substellar companion RV searches. 
There are 62 primaries in total. 

For CBK12, the parent spectroscopic M-dwarf sam¬ 
ple is extracted from a deep magnitude-limited all-sky 
survey, resulting in a large number of M-dwarf targets 
with medium-resolution optical spectra (R ~ 1800, A < 
920 nm). The final subsample of 1452 candidates, deemed 
appropriate for a companion search, is a somewhat 
serendipitous selection from the overall survey (see Sec¬ 
tion 14.1.21) . whose brightness distribution is unimodal 
and peaks at i ~ 16.8. Though the DC13 sample sat¬ 
isfies the colour cuts explicitly imposed by CBK12 (i.e. 
i — z > 0.2, r — i > 0.5), as a whole, the CBK12 sample 
is redder than DC13. CBK12 has median i — z > 0.5, 
mode ~ 0.55 (see Fig lcl therein), compared to median 
i — z ~ 0.4 for DC13, the mode being ~ 0.35. This hints 
that CBK12 samples a somewhat later M-dwarf popula¬ 
tion than our study. 

A known limitation associated with SDSS spectroscopy 
is its fiber collision limit. Two objects within 40" could 
fall within the same fiber, and only one of them would 
have its spectrum taken. If resolved tertiary compan¬ 


ions are commonplace around tight binaries, then only 
a fraction of every such system surveyed would attain a 
spectroscopic measurement of the inner binary and de¬ 
tect an RV variation (Cullen Blake, private communica¬ 
tions). The consequence of missing such systems would 
be an underestimated binarity rate, an effect to which 
CBK12 may be subjected. 

Relative to these former works, the present study 
claims no superiority in the definition of its original sam¬ 
ple selection. As did its predecessors, it makes use of 
an existing dataset whose initial purpose was not neces¬ 
sarily to be conducive to rigorous statistical analysis of 
stellar multiplicity. It is neither volume- nor magnitude- 
limited (see Figure [2] for magnitude distribution of sam¬ 
ple stars). The number of stars monitored by Kepler - 
> 3000 within ~ 200pc- is an improvement over both, 
and samples from both nearby and somewhat more dis¬ 
tant stellar environments., though the corrrection fac¬ 
tors are large for longer periods. Note that the Kepler 
M-dwarf targets are dominated by early-types. 

The differing populations studied by each study (early- 
vs mid- vs late-M-dwarfs) may contribute some intrin¬ 
sic, physical differences between each measure. Af¬ 
ter all, multiplicity is a heavily mass-dependent phe¬ 
nomenon. However, here one expects the difference to 
be small. Even down to the very-low-mass (VLM) and 
substellar regime, the overall multiplicity statistics are 
within a factor of 2 of t heir low-mass star counterparts 
(jDuchene fc Krausil2013ll . In addition, though multiplic¬ 
ity is generally a decreasing function of mass, one would 
expect from scaling and observational evidence that the 
orbital period distribution of companions also moves in¬ 
wards, i.e. the close companion fraction would inflate for 
the VLM stars relative to more massive objects. The two 
effects work in opposing directions. 

4.2.2. Potential Sample Biases 

One concern in surveys with magnitude limits, how¬ 
ever implicit, is that they tend to favour the inclusion of 
unresolved binaries due to their greater combined bright¬ 
ness. In other words, the close binaries are surveyed 
out to a larger volume than the single stars, resulting in 
an overestimate of their numbers relative to a volume- 
limited suryev.^^his pitfall is known as the ‘Branch 
Bias’ rtBranchl I1976D . To counteract this bias, a cor¬ 
rection factor that is dependent on the flux ratio dis¬ 
tribution of the binary components may be applied to 
convert into a volume- limited measurement, w hose an¬ 
alytic form is given in iBurgasser et al.l (120031) (eqn (4) 
and (5) therein). In the case of a strictly magnitude- 
limited imaging survey and assuming all binaries consist 
of equal-mass components, this issue may produce up to 
2.5x upward bias in the bin a rity fraction measur ed (see, 
e.g., IBurgasser et al.l (1200311 : IJddar et~akl (I2013f) l. The 
reality is usually less extreme, though still appreciable. 
For a realistic binary population with overall unresolved 
binary fraction of 20% - 40%, the overestimate inferred 
from the EB sub-population in a magnitude-limited sur¬ 
vey can be ~20 - 60%. Scaling the observed NSPS for 
our period range of interest by this factor reduces the 
best estimate to ~0.07 - 0.09, which is less discrepant 
from the previous measurements. Adding this possible 
systematic bias in quadrature with the formal error on 
the fit, we adopt -0.04 as the lower uncertainty on our 
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NSPS. Note that FM92 make an effort to account for the 
Branch bias by discarding 3 SB systems whose individual 
components do not meet the magnitude selection crite¬ 
rion. CBK12 have not attempted to correct for this bias. 
GAIA parallaxes will be authoritative in the resolution 
of this bias. 

The unresolved nature of these binaries might pose an¬ 
other selection bias for the sample - they may distort 
the colour, hence affect their initial selection and identi¬ 
fication in DC 13, which uses broadband photometry to 
constrain the stellar properties. DC13 in turn makes a 
temperature cut at T e g < 4000AT and a surface grav¬ 
ity cut at log(< 7 ) > 3.6 to restrict their target range to 
low mass dwarf primaries. Bias in the sample can occur 
when M-dwarf binaries tend to be systematically mis- 
classified due to its colour distortion and, as a result, 
are ejected from the sample to begin with. Of course, 
it is also possible that binarity in non-M-dwarf primaries 
may effect a colour change that causes them to be prefer¬ 
entially mistakened as M-dwarfs. We deem both effects 
to be small. On the lower main sequence, luminosity 
scales steeply with mass - as approximately M*. Equal- 
mass binaries should present minimal colour difference 
from the primary itself since the primary and secondary 
would present nearly identical spectral energy distribu¬ 
tions (SED). Unequal-mass binaries should expect the 
primary SED to overwhelm the contribution to the light 
detected overall. 

Though the DC13 stellar parameters are considered 
to be improved upon the original KIC, recent stud¬ 
ies suggest a possible systematic underestimation in 
their radii fe.g. Everett et al .1120lilt iPlavchan et al.lfJOTl : 

iNewton et al.ll2015T) . bv ~15%. If the true radii of the 
M-dwarf primaries are larger than assumed for this anal¬ 
ysis, greater detection probabilities are implied, and the 
weight for each detection decreases by the same factor. 
Such an error could contribute a ~15% overestimate in 
the NSPS. 

4.2.3. Completeness Estimates 

Both former studies have strived to make conserva¬ 
tive completeness estimates. This means their surveys 
are probably more sensitive than assumed for their oc¬ 
currence rate calculations, rendering the results upper 
limits. 

In this work, we have assumed that all eclipsing sys¬ 
tems, even grazing ones, would be easily detected. This is 
a good assumption for stellar eclipses until the extremal 
grazing regime (see Figure ED, hence we may be overop- 
timistic in our detection sensitivity in our formal calcu¬ 
lations. An overestimate in completeness translates into 
an underestimated multiplicity fraction. If this is the 
only factor biasing the completeness, our result would 
actually be a lower bound. 

In the opposite direction, we could be ‘over-complete’. 
That is, perhaps we have mistook a few systems for gen¬ 
uine EBs. It is not far-fetched to expound that our 
‘EB detections’, especially the relatively shallow ones and 
those for which the secondary eclipse are invisible, may 
in reality be due to the transits of low luminosity com¬ 
panions (e.g. brown dwarfs and giant planets). In effect, 
we may worry about the ‘false-positive rate’. 

Existing knowledge of brown dwarf and giant 


planet occurrences as compa nions to M-dwarfs is 


that they are very rare (e.g. Marcv & Benitzj 

1989; 

Endl et al.!i2006: Kraus et al. 2008^ Dieterich et al. 

2012 ; 

Dressing & Charbonneaul 20131; iMorton & Swift 

2014; 

Bowler et al., |2015). According to equation (81 of 

Johnson et al.| (|2010jl which describes the fraction of 


stars hosting giant planets as a function of stellar proper¬ 
ties, the M-dwarf targets deemed to be EBs in our sample 
have less than 4% probability of hosting a giant planet 
within 2.5 AU. Note that one such false positive in the 
original EB candidate list (KID5794240, or KOI-254.01) 
has already been identified, confirmed, and excluded. 
The chance that another giant planet might mimic an 
M-dwarf EB would push the giant planet host fraction 
among these M-dwarfs to 2/14 = 14% within < 0.4 AU, 
which is highly unlikely. RV and imaging campaigns have 
generally found the occurrence of brown dwarf compan¬ 
ions to stars regardle ss of spectral type and orbital sepa¬ 
ration to be ~ 0—5% (iMarcv fe Benitzlll989HKraus et al.l 
l2008t iDieterich et alJl2012HBowler et al.ll2015lf , compris¬ 
ing the so-called ‘brown-dwarf desert’. Among the Kepler 
low-mass stellar targets, only the brown dwarf LHS6343 
( Johnson et al.ll201lt IMontet et al.ll2015l l has been found 
to eclipse an M-dwarf (not included in the DC13 cata¬ 
logue). Still, the likelihood of finding transiting BDs and 
Jupiters in the remaining sample should be very small 
and thus unlikely to affect our results, and future RV 
follow-up of these binary systems will best provide con¬ 
firmation of this assumption. 

A couple of systems in this sample do present a shal¬ 
low primary and lack a secondary eclipse (KID9772531, 
KID11853130). Both have relatively lengthy periods. 
Throughout this analysis, we have been interpreting their 
nature as systems with eccentricity. However, it is also 
possible that one or both of these are false-positives, for 
their eclipse depths suggest companion radii consistent 
with a giant planet or brown dwarf. The rounded bottom 
and depth in the primary eclipse of KID9772531 indicate 
an eclipsing object of radius ~O.O9i ? 0 ~ 0.9i?j, placing it 
on the intersection of giant planet, brown dwarf, and late 
M-dwar f star in the radius-mass relatio n of degenerate 
objects (lBertall2f)l3t IMontet et al .1 1201 Til . KID11853130 
has been included in a most recent compilation of plan¬ 
etary candidates (KOI- 3263.0 1), albei t with a high false 
positive probability dSwift et alil2015ll . 

One system, KID6023859, presents a clear and deep 
( 8 %) primary but the secondary is absent. The fiducial 
period is 27 days but there is some evidence from RV 
data that the actual period could be twice this value 
(Jonathan Swift, private communications). If this is the 
case, the NSPS could be revised upwards by ~0.8%. 

High-resolution spectroscopy and RV monitoring with 
sensitivity to magnitudes r ^16 should be able to elu¬ 
cidate the stellar nature and orbital elements of these 
companions. 

4.2.4. Higher-Order Multiples 

The existence and statistics of an outer hierarchical 
triple companion to short-period binaries have implica¬ 
tions for their formation channels. We do not directly de¬ 
tect transiting higher-order companions among the EBs 
found in this study, hence do not attempt to place a 
constraint on its statistics. Dynamical stability and evo¬ 
lution arguments would suggest that, if triple compan- 
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ions existed, they are likely much farther out with the 
inner binary and not usually coplanar, hence are inac¬ 
cessible via directly observing their transits in a survey 
of the current scope. Nonetheless, other methods may 
be used to infer the existence of further companions. 
For instance, from eclipse timing variations (ETVs), the 
presence of a third, non-transiting companion at ~ 3 
years has been inferred for one of the EBs (KID10979716) 
(jBorkovits et al.f2015f) . KID3830820 is a blend of at least 
3 stars, and these visual neighbours may be physically as¬ 
sociated. The system of KID11853130 also harbours a vi¬ 
sual companion at 0 . 8 ,, J identified through lucky imaging 
(PI Lillo-Box) and adaptive optics (PI Ciardi, see note 11 
to TablcQ]). This EB sample could benefit from more sys¬ 
tematic high-resolution imaging to uncover triple com¬ 
panions at a few AUs. 


5. CONCLUSION 

In the above exposition, we have described an exercise 
in mining data from the precision photometry monitoring 
survey Kepler to conduct a study of M-dwarf multiplcity 
in close (but detached) orbits. The inference is rooted 
in the detection of EBs in this dataset, which has stared 
at ~ 3000 field M-dwarfs as faint as r-magnitude 17 in 
the direction of Cygnus for nearly one year or more. The 
statistical method borrows directly from the calculation 
of exoplanet occurrence rates, with simplifications. Ours 
represents a very simple and straightforward way to in¬ 
fer the likelihood that an apparently single M-dwarf in 
the sky should harbour a companion within a 90-day or¬ 
bital period. It uses a minimalist and transparent com¬ 
pleteness correction, which suffices for its purpose given 
the nature of the technique and quality of the dataset. 
Moreover, this study involves real binary detections on 
the largest M-dwarf sample to date used for this type of 
statistical investigation. 

We find that, among the Kepler M-dwarfs with no re¬ 
solvable companion and whose photometry is consistent 
with that of a single star, there are on average O.lljto '04 
(stellar) companions within a 90-day orbit. Assuming 
all apparently single M-dwarfs are either actually single 
or binary, this translates to a 11 + 4 % binarity fraction. 
This value is significantly higher than previous measure¬ 
ments of M-dwarf multiplicity from literature, though it 
may be partially inflated by potential survey biases (re¬ 
flected in the lower bound error) as well as uncertainty 
regarding the nature of several objects, which could fur¬ 
ther cause 2-4% overestimate. Two studies, FM92 and 
CBK12, both find a value under 4% with 1 — 2% uncer¬ 
tainty for a comparable binary separation regime using 
RV techniques. 

To reconcile the tension between the present work and 
the literature requires careful thought to the comparabil¬ 
ity of the studies, as discussed in Section [I] RV followup 
to the binary detections presented in this paper would 
be critical to confirm or reject the stellar nature of com¬ 
panions and verify the interpretation of these systems. 

Future studies, both planned and in progress, will be 
key to refining this important measurement. One valu¬ 
able venue could come from volume-limited surveys as 
that conducted by the RECONS team, which is a contin¬ 
ual effort to exhaustively map the solar neighbourhood 


(e.g. iHenrv et al.ll2006l h The statistics they will derive 
would form an interesting comparison to the several stud¬ 
ies done on quasi-magnitude-limited surveys, and one 
should look forward to Winters et al. (in prep), which 
aims to quantify the M-dwarf multiplicity issue based 
on the RECONS data. In fact, according to a compre¬ 
hensive private compilation of stars within lOpc of the 
Sun, there is already indication that the early M-dwarf 
population has a close binarity fraction greater than 5% 
( Kevin Apps, priva t e com munications). 

iDeshnande et all (|2013T ) describes an ongoing effort to 
survey more than 1400 M-dwarf stars within the frame¬ 
work of the Sloan Digital Sky Survey, using the near- 
infrared APOGEE spectrograph. Preliminary data prod¬ 
ucts are to be made public as part of the 12th SDSS data 
release. This project is intended to be useful for tackling 
questions of stellar multiplicity, among various other top¬ 
ics. The RV information collected over multiple epochs 
would enable sensitivity to a parameter space in orbital 
period that encapsulates that addressed in the present 
study. 

There is yet another potentially promising route that 
could help settle the issue of close stellar multiplicity, 
using photometric light curve data alone. It involves 
the use of distinct signatures of rotational modulation 
for the detection of multiple systems. An example of a 
possible implementatio n of this concept is described by 
iRappaDort et all (120141 b 

The present study also serves to preview the potence 
and urgency of ongoing and upcoming continuous pho¬ 
tometry monitoring missions li ke t he Kepler 2-Wheeled 
Mission (K2) dHowell et all I2014T) and the Transiting 
Exoplanet S u rvey S atellite (TESS) (|Ricker et al .1120111) . 
iKraus et all (|2015H recently studied in detail a mid- 
M+M eclipsing binary discovered with K2. Undoubt¬ 
edly many will follow. Continuous microlensing surveys 
like the Wide Field I nfrared Survey Telescope (WFIRST) 
(jSpergel et al.ll2013h and t he Korean Microlen sing Tele¬ 
scope Network (KMTNet) (lAtwood et al.ll2012l) will pro¬ 
vide interesting opportunities to study multiple star pop¬ 
ulations in the direction of the galactic bulge. These pro¬ 
grams have a critical role to play in the next decade of 
stellar astronomy. 
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Fig. 7.— The period-folded light curves of M+M eclipsing binary targets that made our final list. Note the dramatic in-eclipse scatter in 
the light curve of KID3830820, which is chiefly due to varying degrees of blended light with neighbours from quarter to quarter. See Table 
|T]for descriptions of each pre-flattened light curve. 
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Fig. 8.— Period-folded M+M EB light curves continued. See Table|j]for descriptions of each pre-flattened light curve. 
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KID 4936334, Period = 0.354369 days 



KID 9077796, Period = 0.23808 days 



KID 12004834, Period = 0.262317 days 



Fig. 9.— Three very tight EB candidates which have duty cycles inconsistent with their theoretical maximum. They exhibit some level 
of contact. KID4936334 has anomalous photometry - its DC13 catalogue entry registers M* = O.O9M0, T e ff = 3433, and Z = —2.0. These 
candidates have been discarded from the final list. All three are Villanova objects. 
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Fig. 10. — Two unconvincing Villanova EBs. They do not meet the criteria for detection by our algorithm, hence are not included in our 
analysis. 














